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Abstract 

We present a pedagogical review of the physics of fractional Chern insulators with a particular focus on the connection 
to the fractional quantum Hall effect. While the latter conventionally arises in semiconductor heterostructures at low 
temperatures and in high magnetic fields, interacting Chern insulators at fractional band filling may host phases with 
the same topological properties, but stabilized at the lattice scale, potentially leading to high-temperature topological 
order. We discuss the construction of topological flat band models, provide a survey of numerical results, and establish 
the connection between the Chern band and the continuum Landau problem. We then briefly summarize various 
aspects of Chern band physics that have no natural continuum analogs, before turning to a discussion of possible 
experimental realizations. We close with a survey of future directions and open problems, as well as a discussion of 
extensions of these ideas to higher dimensions and to other topological phases. 
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1. Introduction 

The application of a transverse magnetic field to a two-dimensional electron gas breaks time-reversal 
symmetry and gives rise to a highly degenerate single-particle spectrum in the absence of electron-electron 
interections: the celebrated Landau levels. When the electronic density is commensurate with the macroscopic 
degeneracy of a Landau level so that an integer number of such levels are occupied, the resulting insulating 
phase exhibits the integer quantum Hall effect (IQHE), with spectacular consequences for charge transport: 
the Hall resistance is quantized to be an integer multiple of h/e 2 , while the longitudinal resistance vanishes 
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[1]. When the density is a rational fraction p/q < 1 of the Landau level degeneracy, there is no appropriate 
single-particle description and it lies to interactions to select a suitable correlated ground state from the 
degenerate manifold. The resulting phases exhibit the fractional quantum Hall effect (FQHE) [2,3]. The 
FQH phases are prototypical example of topologically ordered [4] phases of matter: they are described at 
low energies by a topological quantum field theory written in terms of an emergent gauge held, exhibit 
fractionalization of quantum numbers in the bulk, and have gapless protected chiral edge states. 

The simplicity of the setting in which the QHE was hrst discovered was extremely helpful in unraveling its 
explanation [5,3]. For the two-dimensional electron gases (2DEGs) in semiconductor heterostructures that 
are the venue for all QH experiments thus far, we may ignore most of the typical complications associated 
with the surrounding solid and treat the electrons as having a parabolic dispersion. Lattice effects at best 
lead to a change in the (effective) mass of electrons away from its value in free space. Hand-in-hand with 
this description in terms of nearly free electrons is the fact that the structure of single particle states in the 
lowest Landau level (henceforth abbreviated as LLL) strongly constrains the wave functions of the many- 
body ground states that emerge as a consequence of interactions: they are forced to be analytic functions 
of z = x + iy. The built-in analyticity of the LLL wave functions and the structure of the topological 
quantum held theory that describes their topological order has allowed the full panoply of techniques of 
conformal held theory to be brought to bear in constructing variational wave functions for various QHE 
'filling fractions', v = p/q [3,6,7,8]. Various effective field theory techniques, such as statistical transmutation 
by flux attachment, [9,10] also rely on the simple parabolic dispersion assumed for the electrons. 

Can the QHE be detached from this idealized limit? There are two independent idealizations that must 
be addressed. First, the solid could affect electronic motion more seriously so that the effective mass approx- 
imation gives way to the formation of energy bands; does the QHE survive in this case? Second, a uniform 
magnetic field both breaks time reversal symmetry but also affects electron dynamics at long wavelengths 
in a decidedly unusual fashion: the formation of Landau levels. Are both essential for the QHE? 

The answer to both questions is known for the I(nteger)QHE, which is mostly a single-particle phe- 
nomenon. In a landmark paper in 1982, Thouless, Kohmoto, Nightingale, and den Nijs [11] analyzed the 
uniform-field Hall effect in a strong periodic potential that was known to lead to an intricate spectrum, the 
so-called Hofstadter butterfly; they showed that it gives rise to an integer QHE under certain conditions, 
i.e., whenever the chemical potential lies in a gap. Indeed, the Hall conductance was shown to map to a 
topological invariant associated with filled bands — the (first) Chern number. Six years later, in another 
striking development, Haldane [12] answered the second question, showing by an explicit construction of 
a tight-binding model on a honeycomb lattice that a quantized Hall conductance can arise from a fully 
filled band even in the absence of a net magnetic field. In his model, time-reversal symmetry is broken by a 
spatially inhomogeneous magnetic field with zero average, and the Hall conductance again equals the Chern 
number of the band. 

Of course, the next logical question is: can the FQHE, canonically a property of interacting electrons 
in a fractionally filled Landau level, also be separated from the weak lattice and uniform magnetic field 
limit? The topological flat band models attempt to ask this question in a specific fashion: if it is true for 
independent electrons that a filled Chern band is equivalent to a filled Landau level, then is it also true 
for interacting electrons that a fractionally filled Chern band is equivalent to a fractionally filled Landau 
level? Microscopically, these models proceed to construct a (nearly) degenerate low-energy subspace - the 
'flat band'. While an exact degeneracy requires electron hopping over arbitrarly large distances, the hopping 
amplitudes decay exponentially and thus a relatively flat band can be produced by keeping a small set of 
hopping amplitudes. The suitability of a band for realizing correlated states is naively quantified by its 
flatness parameter: the ratio of the band gap (which sets a bound on the strength of the interactions one 
can safely include) to the bandwidth - although, as we will see, this is only one criterion for a 'good' Chern 
band. 

To understand the significance of a lattice realization of the FQHE it is useful to recall that while in 
theoretical terms the setting in which the FQHE was discovered is quite simple, experiments must go to 
great lengths to achieve a limit where this simplicity emerges. First, a variety of growth techniques must be 
employed to confine electrons in a two-dimensional plane, while simultaneously maintaining a high electron 
mobility - the latter restriction because too much disorder destroys the FQHE. Second, an extremely strong 
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magnetic field - of the order of tens of tesla - must be applied for Landau quantization to be appreciable. 
Finally, the systems must be cooled to below at least 10 K to reach the energy scales at which FQH physics 
is manifest. All in all, this involves a tour de force in experimental technique, and reliably reproducing the 
full sequence of FQH states remains a significant challenge even three decades after the original discovery. 
In contrast, in a lattice model, the characteristic energy scales can be in the tens or even hundreds of kelvin 
(assuming a lattice energy of a few meV) and are naturally in the limit where each unit cell sees a large 
effective magnetic field. Therefore such systems may realize robust phases that are concomitantly insensitive 
to disorder and thermal fluctuations. The possibility of stabilizing exotic topological phases outside of a 
dilution fridge and without a superconducting magnet, and the avenues it opens for further experimental 
probes of these phases, is among the primary motivations of the activity in this field. On a more theoretical 
front, the interplay of lattice symmetries with topological order leads to new physics, unique to the lattice 
realizations of the FQHE. 

The recent flurry of interest in the FQHE in topological flat bands began with the work of three groups 
[13,14,15]; several authors have since constructed models with nearly flat (non-dispersing) Chern bands 
which allow interactions to dominate at partial fillings while leaving the gap to neighboring bands open. 
Initial papers [15,16,17] reported evidence for FQH states at v = 1/3, 1/5 and (for bosons) 1/2, in finite-size 
studies of short-ranged interactions projected to these bands. A detailed study of properties of ground state 
entanglement [18] allows a rationalization of some aspects of these results in terms of a generalized Pauli 
principle familiar from the LLL problem. The connection to the FQHE was made more explicit by Qi [19] 
who gave a fairly general recipe for translating familiar model wave functions and Hamiltonians from the 
LLL to Chern bands on cylinders by elegantly mapping Landau gauge eigenfunctions to particular Wannier 
functions. Finally, work by the present authors [20,21] gave a mapping between Hamiltonians in the LLL and 
those in the Chern band by studying the correspondence of the operator algebras that arise in the two cases. 
The outcome of these three early approaches is to identify three desiderata for a Chern band: (i) that it be 
energetically flat with a large gap and thus a large flatness parameter, so that the low-energy description 
may be approximated by projecting into the band (ii) that it be Berry flat, i.e. have near-uniform Berry 
curvature, so that the algebra of densities projected to the Chern band resembles that of similar operators 
in the LLL and (iii) that it satisfy a set of conditions on the Fubini-Study metric of band eigenstates that 
quantify the connection to the LLL. Building on this, a reformulation of the Hamiltonian theory of the 
FQHE has been applied to study Chern bands [22,23]. 

This is as good a place as any to note that on a lattice the distinction between having a net magnetic field 
and not having it at all is not as sharp as it may seem. Essentially, it is always possible to stick a full flux 
quantum through some subset of loops on the lattice to shift the average magnetic field without affecting 
the actual physics. From this perspective, the physics in these flat-band models has a family resemblance to 
earlier studies of lattice versions of the FQHE [24,25] with uniform magnetic fields. In this earlier work, the 
authors studied a fixed filling factor while varying the flux per plaquette from small values and large unit 
cells, where the standard Landau level description holds, to somewhat larger flux values and smaller unit 
cells, where that description broke down. As they were able to change this parameter without any evidence 
of encountering a phase transition, the latter limit constituted an observation of the FQHE in the presence of 
strong lattice effects. As we demonstrate below, a more analytic approach clarifies the equivalence between 
this earlier Hofstadter and the current Haldane versions of the FQHE. Furthermore, this connection permits 
the construction of closed algebra of one-body operators [20,21,22,23] that gives a Chern band realization 
of the Girvin-Macdonald-Platzmann algebra [26] of the LLL. 

In addition, several authors have examined aspects of Chern band physics where the presence of the lattice 
introduces an additional layer of richness and complexity. Wc discuss two examples in this review. The first 
of these is the intriguing question of the transition from a kinetic-energy dominated regime where the inter- 
particle interactions are weak and the system is in a gapless metallic phase (superfluid, if bosonic) to the 
strong-coupling regime where the flat band FQHE sets in. Continuous bandwidth-tuned transitions between 
these two limits provide instances of topological phase transitions [27,28], which we discuss briefly below. 
A second example of the reappearance of lattice physics is the case of bands with higher Chern number. 
Closely related to 'bilayer' quantum Hall states in a Landau level, higher Chern bands both provide a simple 
route to non-abelian FQH states [29] . In addition, the lattice physics reappears via a nontrivial interplay of 



3 



crystalline point defects with the topological order, endowing them with a nontrivial 'quantum dimension.' 

At the time of writing, FQH physics in Chern bands is numerically well-tested, with the proposed can- 
didates having passed a host of detailed checks by various groups. However, an experimental realization of 
a fractional Chern band remains elusive, although several candidates have been suggested. We will discuss 
a few promising ones below: these include solid-state realizations in oxide interfaces [30], as well as optical 
lattice examples both with short-range [31] as well as dipolar interactions [32,33]. 

1.1. This Review 

Before proceeding, we list a few disclaimers as to the scope of this review. As with any rapidly progressing 
field, there is a risk that some of what we discuss will be out of date by the time of publication, or that 
currently active directions fail to realize their early promise. To avoid this problem of obsolescence and in 
the interests of pedagogy our primary focus is on the case where the characterization at the present time 
is most complete: Chern bands in two dimensions. We will eschew a discussion of other topological bands 
- such as those relevant to fractionalized analogs of time reversal invariant Z 2 topological insulators - 
and higher dimensions. The bulk of the review will focus on bands with Chern number C = 1, although 
we will briefly discuss bands with C = 2. We will also restrict our detailed discussion to the early work on 
the subject, and devote a significant fraction of this review to elucidating the relationship between FQHE 
states realized in the continuum, and their putative analogs in Chern bands. We will attempt to present a 
broad summary of the subsequent direction of the field, with pointers to the literature for the interested 
reader. We will also briefly discuss new aspects of FQHE physics unique to Chern bands, as well as potential 
experimental realizations. Also, while we describe the important numerical results, we will not discuss the 
specifics of the numerics in this review, referring the reader to the original literature for further details. Our 
intention is that by the end of this article, a novice to the subject will emerge familiar with the broader 
aspects of fractionally filled topological bands, and equipped to tackle the (rapidly growing) literature; we 
trust that the reader will not find our approach too idiosyncratic. 

2. Survey of Models 

We begin with a brief introduction to Chern insulators in general and flat band models in particular, 
that will serve both to fix notation and provide background for the remainder of this review. Following a 
general discussion of tight-binding models without interactions, we proceed to the question of interacting flat 
band models where we review the numerical results that initiated the systematic study of fractional Chern 
insulators. We also briefly discuss measures such as ground state entanglement, that can be unambiguously 
used to identify the topologically ordered phases in these models. 

2.1. Lattice Models of Chern Insulators 

Consider a tight-binding lattice model of the form 



where the {t^} are the (in general complex) hopping matrix elements, the i,j label Bravais lattice sites, 
and a,b = 1, 2, . . . , Af are internal indices that label different orbitals or sites within a unit cell. Invoking 
Bloch's theorem and working in momentum space, we may write the Hamiltonian as 



where a, b are indices that label sites within a unit cell, and k is the crystal momentum restricted to the 
first Brillouin zone (BZ). (Here and below, we will explicitly indicate when repeated indices are summed 
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over, and unless otherwise specified will work in two dimensions.) The solution of the Af x M eigenvalue 
problem J^b ^ab(k)u^(k) = e Q (k)u"(k) defines the Bloch bands e Q (k) and Bloch states u"(k). We will take 
the corresponding eigenvectors to be normalized, ^ a |u"(k)| 2 = 1. The corresponding eigenstates are given 

by 

|k,a)= 7 Ul°) = E<( k ) c Ll ) ( 3 ) 

a 

and in terms of the operators 7^ a we have H a — J2^=i e a(k)7k Q 7k,a- 

The Chern number of a given band a is a topological invariant which can be defined only if the band is 
isolated from all other bands, and is computed as 



X - J^d 2 kB a (k). (4) 



Here, B a (k) is the Chern flux density (Berry curvature), defined as the curl of the Berry connection (Berry 
gauge potential), B a (k) = Vk x «4. Q (k). In terms of the Bloch states, we have 

N 

•4a(k) = *5>r(k)VkU?(k). (5) 

6=1 

A filled band with Chern number C a yields a Hall conductance <jh = C a e 2 /h regardless of whether it arises 
in a system with a net magnetic field [11] ("Hofstadter band") or zero net magnetic field [12] ("Haldanc 
band"). We shall refer to both as Chern bands. In addition, we will assume that we are considering bands 
with C a = l unless otherwise specified. 



2.2. Engineering Flat Bands 



A Chern band shares one important feature with a Landau level: a nonzero Chern number. A Landau 
level has the additional feature that it is an exactly degenerate manifold, so that the only energy scale 
within a fractionally filled Landau level is provided by the interparticle interactions. In a typical Chern 
band, however, the kinetic energy from the band dispersion is a competing energy scale. In order to increase 
the efficacy of interactions the kinetic energy can be 'quenched' by flattening the dispersion, as was noted 
by several authors; it is important that the single-particle Hamiltonian remain local even after the band- 
flattening procedure. Consider a lattice model of the form above, and let us assume the lowest band 1 has 
Chern number 1 and gap A = min(e 2 (k)) — max(ei(k)). Simply performing the transformation 

Mk)-»Mk) = Mk)/ Cl (k) (6) 

and Fourier transforming to return to real space defines a tight-binding model in which the hopping matrix 
elements decay asymptotically as if!- oc e - !* - - 7 '!^ and the lowest band is flat. Thus, as a point of principle, any 
energy band that is gapped away from all other bands can be made flat, at the cost of adding exponentially 
decaying longer-range hoppings; while the resulting tight-binding model is not strictly local (in the sense of 
having matrix elements with compact support), most arguments based on topological order remain valid for 
the exponentially local Hamiltonian that results in this fashion. Of course, truncating the hopping elements 
for sufficiently large \i — j\ renders the Hamiltonian strictly local at the cost of introducing some dispersion, 
characterized by a bandwidth t. Correlations then play the dominant role only if the interaction scale is 
much larger than the bandwidth, U >• t. The maximum interaction scale is itself limited by the band gap, 
since the approximation of restricting attention to the lowest band breaks down unless [/« A. Therefore 
an appropriate quantitative measure of the suitability of a Chern band for realizing correlated phases isthe 
flatness parameter / = A/t. Particles in bands with large / have correlation energies much higher than 

1. The assumption of the lowest band is for pedagogical reasons; the generalization to the case when there are several filled 
bands below the one of interest is trivial and left as an exercise to the reader. 
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their kinetic energies, for interacting Hamiltonians that may safely be approximated by their projection 
into the band. Several authors [34,35,36,37] all constructed flat bands with fairly large flatness parameters, 
even with hopping restricted to a few near neighbors. This solves the first problem we noted earlier: we 
have successfully produced examples in which it is reasonable to approximate the physics by projecting to 
a single band with non-trivial topology. 

2.3. Interacting Flat Band Models 

Once a noninteracting Chern insulator hamiltonian H$ has been obtained that hosts a flat band a with 
large flatness parameter, the next step is to add inter-particle interactions. The full interacting hamiltonian 
takes the form H = H a + V where V is an interaction term, which typically (but not always) takes a 
generalized Hubbard form: 

V=^^2u ij n i n j (7) 

where hi = c\a is the number operator for electrons on site i. For such a density-density interaction, assuming 
a translationally invariant Hamiltonian, we have C/,j = E7(rj — Vj) = J (2?r) 2 U(q)e' lcl '( Vi ~ r i\ The interaction 
term can then be written V — ?7(g)n q n_ q where n q is the Fourier transform of hi. As a measure of the 
interaction scale U we can take (for instance) the maximum value of U (r 4 — r^) . Note that the full space of 
interacting Hamiltonians is somewhat more general than the density-density interaction assumed above - 
for instance it could include ring exchange, pair hopping and other non-density contributions. Nevertheless 
we will mostly restrict ourselves to such Hubbard-type models as they are a reasonable starting point and 
have a well-founded physical origin in an onsite Coulomb repulsion. 

2.4. Projecting to the Flat Band 

Since we assume U -C A we may safely neglect the mixing between a and the remaining bands. This 
allows us to project to the partially filled band; such a projection is implemented by the operator V a = 
|k, a)(k, a\. In second-quantized form, the resulting Hamiltonian takes the form 

H = V a HV a = e «(q)7 q ,«7q, Q + \J2 U (q)Pq;«P-q;a (8) 
q q 

where we have defined the projected density operator in the Chern band, 

Pq;a=PaW = £ U ?*( k+ f) U ?( k " f) >< 7^ , a 7k- f ,a (9) 

If the bandwidth is small compared to the scale of the interactions, t <C U, e a (k) may be treated as 
constant and thus ignored. With this approximation, we finally arrive at the low energy effective flat band 
Hamiltonian, which takes the form 

^cl.^E^W-q;* (10) 
q 

One encounters a similar Hamiltonian in the treatment of interactions in the LLL in a large magnetic field. 
In that case, the effective Hamiltonian of a clean system obtained by projecting density-density interactions 
to the LLL has the form 

Hlt^l^V^e-^Xp^ (If) 
q 



G 



where in this case p q differs from the projected density, Vpq'P by a g-dependent constant, VpqT 1 = 
e -9 2 ^fj/ 4 e j q- R = e~ q2e ^^ 2 p (i where V is the operator that projects to the LLL, R is the 'guiding center' 
coordinate and £b = (Hc/eB) is the magnetic length. V(q) is once again the Fourier transform of a two- 
body density-density interaction. 

The similarity between Hamiltonians projected to a Landau level and to a Chern band, coupled with the 
fact that both the LL and the Chern band have a nontrivial Chern number suggests that similar fractionalizcd 
phases may arise in the Chern band as in the FQHE. The missing piece in making this analogy precise is 
to match the commutation relations of the projected densities, which we discuss at length below. First, we 
give a brief overview of the numerical studies that have driven much of the early progress in the field. 

3. Numerical Results 

Once a flat band Hamiltonian has been constructed, the next logical step is to ascertain the properties of 
its ground state at a given filling. Analogously to the study of the FQHE in the LLL, the approach of choice 
is to exactly diagonalizc the flat band Hamiltonian for small numbers of electrons, and try and extract from 
this the ground state properties in the thermodynamic limit. Since such analysis is built from the existing 
intuition gleaned from LLL numerics, it is useful to first review those original aspects. We should preface 
this discussion by noting that while much of the numerical intuition for the LLL is obtained in the spherical 
geometry where even topologically ordered ground states are non-degenerate, here we discuss continuum 
results that study FQH states on a torus as they can be directly compared to lattice Chern bands. 

On a torus, FQH states are expected to exhibit a ground state degeneracy in the thermodynamic limit: 
they should have a low-lying multiplet of topologically degenerate states, e.g. there are m such levels for 
the 1/m Laughlin state. 2 The gap between these and the remaining states in the spectrum should scale 
to a constant in the thermodynamic limit. Furthermore, upon 'flux insertion' through a handle of the torus 
- implemented operationally by diagonalizing the problem with twisted boundary conditions - the states 
should exhibit 'spectral flow' within the low-lying multiplet: as the flux is increased, the levels should evolve 
into each other, with a periodicity in m flux quanta for the 1/m Laughlin states. For each degenerate ground 
state, a many-body Hall conductance can also be computed by suitably twisting boundary conditions, and 
is expected to match that of the FQH phase. An additional test is to study the uniformity of the particle 
densities, as the FQH state is expected to be a uniform incompressible liquid. The overlap of a numerical 
ground state with trial wave functions — which can be analytically demonstrated to be in a particular FQH 
phase - provides further characterization of its properties. Similarly, several of the FQH wave functions are 
gapped exact ground states of local model 'pseudopotcntial Hamiltonians'; adiabatic continuity between 
these and the actual Hamiltonian can be used to establish the ground state topological order. The fractional 
statistics of the quasiparticle excitations of the FQH also follow a 'generalized Pauli principle' that governs 
the structure of low-lying many-body eigenstates with a fixed number of excitations in a precise manner 
that can be tested numerically. Finally, a more sophisticated tool is to study ground state entanglement, 
whose properties can be rationalized in terms of topological field theoretic considerations and the bulk-edge 
correspondence of the FQHE. Each of these approaches, to a greater or lesser degree, has been transcribed 
to the problem of fractionally filled Chern bands. 

3.1. First Results: Overlaps, Hall Conductances and Spectral Flow 

The earliest evidence for FQH physics took the form of studying ground state degeneracies as well as 
spectral flow, as well as a computation of ground state Hall conductance 3 For instance, Neupert et. al. [15] 



2. However, it is possible that other states, with broken symmetries - for instance those with charge density wave order- 
also exhibit a similar degeneracy on a torus as first noted by Haldane [38] for the case of continuum LL problem, and so some 
care must be taken in using the ground state degeneracy alone to judge if the ground state is in a FQH phase. 

3. We stress that this was at that point state-of-the-art, since the analytic understanding of the FQH in Chern bands was 
then insufficient to implement various other checks. 
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Figure 1. Results from Ref. [15], illustrating (a) energy spectrum for 6 electrons at v = 1/3; (b) spectral flow within the 
ground state multiplet in the fractional Chern insulator; and (c) the lack of such spectral flow in a fractionally filled trivial 
insulator. 



studied systems with 6-electrons at 1/3 filling, and demonstrated spectral flow between the three lowest 
eigenstates and an absence of flow among higher levels. Although the gap between some of the putative 
ground states and the remainder of the spectrum was not particularly large and comparable to the splitting 
between the three lowest states, the spectral flow coupled with numerical estimates of the Hall conductance 
provided early evidence for topological order. Subsequent numerics by Sheng et. al. [16] - on a slightly 
different flat band model, slightly larger system sizes, and appropriately tuned interactions - found gaps at 
fillings 1/3 and 1/5 which were much larger than the spread of the lowest lying multiplet 2. Furthermore, 
while exact parent Hamiltonians in the Chern band were still unknown, the Hamiltonians realizing these FQH 
phases were of the same general form (in terms of range of repulsion, etc.) as the Haldane pseudopotentials 
in the LLL - for instance, to stabilize the 1/5 state, next nearest neighbor interactions needed to be added. 
This is perhaps a good place to comment on the use of Hall conductance as a probe of topological order: 
as demonstrated first by Kol and Read [39], in the absence of Galilean invariance - for example, due to the 
presence of a lattice potential - the connection between the filling and the Hall conductance is somewhat 
more subtle, particularly for more complicated fractions, such as 2/5. 4 



3.2. Quasihole Counting and Generalized Exclusion Statistics 

A careful study of the quasihole excitations of fractional Chern insulators was used by Bernevig and 
Regnault [18] to argue extremely convincingly for the existence of FQH states in the Chern band. After 
first establishing that the gap, ground state degeneracy and Hall conductance all remained robust with 
increasing system size, and that the ground state wave functions described uniform-density fluids (as opposed 
to modulated charge density waves), they proceeded to study the excitations of the putative FQH ground 
state. In additional to being fractionally charged, these are expected to exhibit fractional statistics both in 
their exchange with each other, as well as in the generalized 'exclusion principle' that governs the counting 
of their low-energy many-body eigenstates [40]. While fractional charge and exchange statistics are difficult 
to study in finite-size numerics, the counting of quasihole states can be directly verified. 

Assuming the system at commensurate filling to be in a particular FQH phase places strong constraints 
on the precise number of low-lying energy eigenvalues of the quasihole states in each generalized momentum 
sector. In the conventional FQHE, this counting can be rationalized in terms of a generalized exclusion 



4. This was revisited in the context of fractional Chern insulators by Shankar and Murthy [22,23]. 
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Figure 2. Low energy spectrum of fractional Chern insulators and other phases in fractional Chern bands from 
Ref. [16]. The parameters U and V control the scale of the nearest neighbor and next-nearest neighbor repulsion terms in the 
Hamiltonian. The top panels in (a) and (b) clearly illustrate the difference in scale between the energy spread of the ground 
state multiplet and the energy gap. 

principle [40]. For example, working in the one-dimensional basis characteristic of the continuum Landau 
levels, the generalized exclusion principle for the 1/m Laughlin states forbids more than 1 particle in m 
adjacent orbitals. This leads to a precise prediction for the counting in each momentum sector (here by 
'momentum', we mean the appropriate conserved quantum number conserved by the symmetries of the 
chosen gauge-fixing.) Such a principle can be applied to each conserved momentum sector, and at a more 
coarse grained level the total number of low-energy states in a finite system with a given, fixed number of 
quasiholes is also given by the exclusion principle. Note that these statements are true for the full spectrum 
only in the case of model Hamiltonians designed to render a particular FQH wave function and its quasiholes 
exact eigenstates. For a generic Hamiltonian, they are expected to apply to the low-lying energy spectrum 
that is believed to capture the universal topological content of the phase, which is typically separated by a 
gap from spurious content at higher energies - although such a gap may be small or absent for small system 
sizes or infelicitously chosen aspect ratios. 

In order to rationalize the results for the Chern band, Bernevig and Regnault [18,41] gave an 'unfold- 
ing' prescription to obtain a one-dimensional generalized momentum parameter to label the many-body 
eigenstates: for an N x x N y system, for instance, they labeled states in terms of k = k x + N x k y . With 
this procedure, it is possible to compare the counting of quasihole states both in total as well as within 
each momentum sector; such states were generated by diagonalizing the problem at a filling which deviates 
slightly from commensuration 5 , corresponding to studying the system with a fixed number of quasiholes: 
for instance, diagonalizing the problem of 9 electrons on a 5 x 6 system introduces 3 quasiholes of the 1/3 
state. The total counting matched that predicted for the LLL quasihole states. Although the counting by 
momentum sector is somewhat more sensitive and occasionally exhibited a mismatch, in many cases this 
also conformed to the expectations from the LLL. This strongly bolstered the original numerical studies and 



5. The small system size means that this choice can be delicate; see [18] for details. 
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Figure 3. Results from Refs. [18,41] for N = 8 bosons on an N x = N y = 4 lattice (y = 1/2 bosons). Clockwise from 
top left: (i) kagome lattice model, with tunable hopping phase ip; (ii) (a) energy gap (A^) and finite-size splitting of degenerate 
eigenstates (5) and (b) entanglement gap (At), as function of ip; (iii) particle entanglement spectrum keeping Na = 4 bosons, 
showing entanglement gap (counting of levels below gap match predictions from quasihole counting); and finally (iv) low-lying 
energy spectrum at v = 1/2 showing ground state degeneracy (magnified, inset) and spectral gap. 



provided the first evidence that the FCI excitations obeyed similar generalized exclusion statistics as their 
cousins in the LLL. 



3.3. Entanglement Spectrum 

A powerful tool that has emerged in studying correlated topological phases, particularly in numerics, is the 
study of entanglement properties of the ground state [42,43]. Various aspects of ground state entanglement 
have been discussed in the literature, and it is well beyond the scope of this review to give an exhaustive 
summary; here we will only provide a telegraphic account of the details needed to discuss the numerics. The 
study of the entanglement encoded in a ground state |\E') proceeds first by constructing the density matrix 
of the state, p = |^)(^|. The next step is to divide ('cut') the Hilbert space into two disjoint subsystems A 
and B: Sj = S^a ® F)b, and construct the reduced density matrix of one subsystem, say A, by tracing over 
the degrees of freedom in the remainder: 



p A = Tr s [p] 
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(12) 



The reduced density matrix captures the various properties of ground state entanglement. If we define 6 
an 'entanglement Hamiltonian' Ha we may write pA = e~ HA , and the set of eigenvalues of Ha constitute 
the entanglement spectrum of the ground state [44]. Another commonly used measure is the entanglement 

entropy [45], which is the von Neumann (S v n — Tta[pa hi Pa]) or Renyi (S^ = Tr^p^]) entropy computed 
from the reduced density matrix. 

Devoid of any other information, the entanglement spectrum does not offer significantly more information 
than the entanglement entropy; typically the best that can be hoped for is to understand the level statistics 
of the eigenvalue distribution. Where a knowledge of the entanglement spectrum comes into its own is 
when the system is endowed with an additional symmetry that can serve as an organizing principle for the 
entanglement eigenvalues; if the 'cut' commutes with the symmetry, the entanglement spectrum is block 
diagonal in the corresponding symmetry sectors, and various additional properties can be extracted. For 
instance, in pioneering work Li and Haldane demonstrated that the entanglement spectrum of rotationally 
invariant quantum Hall states on the sphere could be organized using their conserved angular momentum. 
They matched the counting of the low-lying entanglement levels 7 in each angular momentum sector with 
that of the chiral conformal field theory that describes gapless excitations at the physical edge of the 
corresponding FQH state. An enormous body of knowledge about the entanglement spectrum of FQH states 
in the LLL has been built by analytical and numerical studies. 

Bernevig and Regnault [18] computed the ground-state particle entanglement spectrum - which is obtained 
my making a 'cut' in particle number space, i.e. producing a reduced density matrix by tracing out a 
certain number of particles rather than over a geometric region. 8 The resulting entanglement spectrum 
was organized using a similar 'unfolding' procedure as that used for the true spectrum. While the real- 
space entanglement spectrum is tied to the behavior of edge excitations of a FQH droplet, a particle space 
cut encodes the structure of quasihole states (at a filling corresponding to the number of particles left 
after the cut is made.) Thus, the particle entanglement spectrum should exhibit similar counting as the 
actual quasihole spectrum for the same number of quasiholes (see Fig. 3, bottom right, for an example). 
Verification of this property, which has been demonstrated in detail in the LLL, for the fractionally filled 
Chern band constitutes another important indication that the two systems exhibit the same topological 
order. Furthermore, both the energy and entanglement spectra can be rationalized more precisely with an 
understanding of the emergent many-body symmetries (originally introduced by Haldane for the LLL) in 
the context of the Chern band [41]. 

3.4. Extensions: Trial wave functions, Pseudopotentials, and Adiabatic Continuity 

Since the early results that established the existence of FQH phases in Chern bands, there has been 
an explosion of numerical studies on fractional Chern insulators, an exhaustive list of which would be 
impossible to provide. Limitations of space prevent us from discussing several interesting extensions, such as 
to bosonic systems, the inclusion of disorder and numerical studies of bands of higher Chern number. 9 One 
development we do wish to briefly summarize, especially as it connects naturally to the next section, is the 
study of trial wave functions and pseudopotential Hamiltonians. 10 Here, Qi provided a route to constructing 
pseudopotential Hamiltonians and corresponding model wave functions in the Chern band that built on the 
understanding of those in the LLL via a mapping of Wannier states, discussed below. Subsequent detailed 



6. Note that p\ is positive definite so Ha is Hermitian. 

7. Here and wherever we discuss entanglement spectra, by 'low-lying' we mean states below the so-called 'entanglement gap' 
that is generically observed to separate the universal content of the entanglement spectrum from spurious behavior at higher 
entanglement 'energies'. 

8. Note that in a periodic system with topological order there is some ambiguity as to the treatment of degenerate eigenstates; 
Bernevig and Regnault chose to study a uniform incoherent superposition of the q degenerate Laughlin states, but the phase 
dependence of the entanglement entropy of coherent superpositions can be used to extract additional topological properties [46]. 

9. Although we briefly discuss a variational wavefunction analysis for phases in C = 2 bands, below, and comment on some 
analytical results on lattice defects in such bands. 

10. We present this discussion here for reasons of organization; to the uninitiated, its clarity may be much improved by 
returning to it after reading the next section. 
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work extended this pseudopotential formalism [47] , studying the effects of inhomogeneous Berry curvature, 
as well as providing various prescriptions for improving the trial wave functions to account for the detailed 
structure of a Chern band. Finally, adiabatic continuity between Qi's model Hamiltonians and more realistic 
projected density-density interactions in a Chern band has been established numerically for certain FQH 
phases [48]. 



4. Landau Levels and Chern Bands: A Critical Comparison 



As noted in the previous section, a variety of different models realize topological flat bands that host 
insulating phases at fractional fillings. On one level, this is extremely satisfying as it vindicates a certain 
prejudice: that since the IQHE has a lattice realization in a Chern band, so should the FQHE by analogy. 
This is bolstered by evidence such as the fractional Hall conductance, that are consistent with FQH physics. 
However the evidence we have reviewed so far has been predominantly numerical; it is clearly desirable to 
have a more analytical approach underpinning our understanding of the connection between Landau levels 
and Chern bands. Ideally this would achieve two goals: 

i establish some notion of adiabatic continuity between the FCI ground state and model FQH wave 
functions in a Landau level. 

ii provide a quantitative measure of the role of the lattice potential in the FCI. 

To see why such an analytical understanding is important, note that the restriction of the problem to the 
low-energy Hilbert space corresponding to the lowest Chern band follows from the first of our three criteria 
for a 'good' Chern band. In simple terms, we have shown that it is reasonable to study the problem restricted 
to a space of states that have the same global topology (as quantified by their Chern number) as the orbitals 
in a single Landau level. If we can also demonstrate that a simple local Hamiltonian takes a similar form 
when written in the basis of Chern band eigenstates as in the Landau orbital basis, then we can argue 
using the basic notion of adiabatic continuity that the two problems represent essentially the same physics. 
There are two ways to accomplish this: the first, taken by Qi, is to construct a one-one-mapping between 
Landau gauge orbitals and (linear combinations of) Chern band eigenstates, and transcribe a Hamiltonian 
appropriate to a particular FQH phase from the LLL to the lattice Chern band. If this mapping is local, this 
would establish the existence of an FQH phase in the Chern band. An alternative approach, put forth by the 
present authors, is to begin with the same model of a projected density-density interaction in both cases, 
and demonstrate that the algebra of projected density operators in the Chern band has the same features 
as the Girvin-Macdonald-Platzman algebra in the LLL. While the former approach permits a construction 
of trial wave functions, an advantage of the latter approach is that makes manifest the importance of the 
second criterion of Berry flatness for a Chern band to reproduce LLL physics. A more detailed study of the 
algebra also suggests the various constraints on the structure of the single-particle states in the Chern band 
that constitute the third and final criterion. 



4.1. Wannier Orbitals, wave functions, and Pseudopotential Hamiltonians 



The first concrete correspondence between Chern bands and Landau levels was proposed by Qi [19], who 
constructed lattice analogs of the Landau gauge single-particle eigenstates familiar from the continuum. 
Recall that the latter are given by 

1M') = 1/4 ^= e^" e -i%(*-^)' (13 ) 



which are centered at x = K y £ B , and K y = with n an integer. Clearly, the center of this wavefunction 
shifts as x — > x + j^i 2 B under transforming K y — >• K y + ^ . In the Chern band, Qi constructed Wannier 
states that have a one-dimensional momentum label with respect to which their centers evolve with a similar 
'shift' property We briefly summarize his Wannier construction for a Chern insulator placed on an infinite 
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cylinder, periodic in y and infinite in x. We will follow the notation of the previous section and consider a 
Chern band a. It is always possible to perform a unitary (gauge) transformation on the single-particle states 
so that we fix one component of the Berry vector potential: (A a ) y — 0. In this case the Wannier states in 
the Chern band are given by (suppressing the index a) 

\W(ky,x)) = ±- y^e-iJo^^^e-i^-Oiky)/^)^) (14) 
Lx K 

where k = (k x , k y ), 0(k y ) = J Q 27r A x (p x , k y )dp x , and x is an integer labeling the lattice sites. The phase factor 
involving 6(k y ) guarantees that the symmetry of the Bloch functions under k x — > k x + 2tt is properly ac- 
counted for. Under a gauge transformation that multiplies the Bloch states by a phase factor |k) — > e 4< ^ k )|k), 
the Wannier function is invariant upto a phase: \W a (k y ,x)) — > e l ^°' kv ^\W(k y ,x)). Most significantly, the 
center-of-mass position of the Wannier function \W(k y ,x)) is given by 

(x) = (W(ky,x)\x\W(k y ,x)) = x- 9(k y )/2w (15) 

in other words, 6(k y )/2n is the shift of the Wannier function from the lattice site, often referred to as the 
charge polarization. However, the integral of this quantity gives the Chern number, C = — ^- L 71 6{k y )dk y ; 
thus, if C 0, the Wannier function is not periodic in k y . The center-of-mass position shifts from x — > x + C 
when k y evolves from to 27r, so that the Wannier functions satisfy twisted boundary conditions: 

\W{k y + 2ir,x)) = \W(k v ,x + C)) (16) 

We can define a generalized momentum coordinate K y — k y + 2Trx for < k y < 2ir, in which case if we define 
\W{K y )) = \W(k y ,x)), we have a continuously defined Wannier function for K y e R. The center-of-mass 
position of the Wannier function now evolves continuously for K y , which thus plays a role analogous to the 
conserved momentum of the Landau gauge eigenfunctions in the continuum. It may be verified that under 
the weak-field (Hofstadter) limit Ib ^> 1, the Wannier functions (r\W(K y )) — > ip Ky (r). 

The next step is to recall that many model wave functions representing quantum Hall states in the LLL are 
exact ground states of local pseudopotcntial Hamiltonians [6,49]. Working on the cylinder, we can represent 
the model wave functions in the basis of occupation numbers of the single-particle states (13); using the 
correspondence demonstrated above, this can be transcribed directly into a wavefunction on the lattice 
written in the basis of Wannier orbital occupation numbers. For instance, consider the LLL wavefunction 
for the m th Laughlin state on the cylinder 

^({z i }) = ({z i }\^) =Ql[(e^ L * _e 2 ^/^) m e-E-?/^ (17) 

i<j 

where f2 is a normalization factor and z — x + iy. Following Qi's prescription, the corresponding Chern band 
Laughlin wavefunction on the lattice takes the form 

l*c / B m ) = E $ (K})III W W^) (18) 

where the coefficients are obtained by expanding the wavefunction in the basis of Landau gauge eigenstates: 

^Ki) = ^ I n^^/L y W<T({^})- (19) 

The pscudopotential Hamiltonian that renders (17) an exact ground state can also be written in terms of the 
single-particle orbitals, and therefore readily transcribed to the Chern band. A final question is to determine 
if the resulting interacting lattice model is local. It is possible to demonstrate that the model Hamiltonian 
is indeed built from local terms on the underlying lattice. In this fashion, Qi produced a model Chern 
band Hamiltonian and corresponding exact ground state that has the same topological order as a given 
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quantum Hall wavefunction, thereby establishing the existence of the FQHE in a Chern band. In essence, 
the pseudopotential Hamiltonian is a special point in the space of Chern band Hamiltonians where the FQHE 
is manifest. It is not clear a priori if this is related to the projected density-density interactions typically 
used in lattice models. The trial wave functions that result from this construction may be compared with 
the exact ground states of projected Hamiltonians; in doing so, some care is needed to gauge-fix the single- 
particle eigenstates in a manner that preserves the underlying lattice symmetries. For a detailed discussion, 
we refer the reader to the literature [50]. 



4.2. Projected Densities and Woo Algebras 



An alternative perspective to Qi's Wannier orbital approach is to consider the algebra of densities pro- 
jected into the Chern band which will allow a quantitative comparison between the LLL and the Chern 
band [20]. Recall that in the quantum Hall problem, the magnetic translation operators 11 obey the so- 
called Girvin-Macdonald-Platzmann Algebra [26], 

K 1) P q J=2-in(^^l)p qi+q2 (20) 

Consider the projected densities in a Chern band, given by (9); at long wavelengths qa <C 1, we may 

■ f k+q/2 dV A 

expand (k + § ) v% (k - § ) « 1 - iq ■ E6«6*( k )^ u 6 ( k ) ~ e J "-"/' ' a{ so that (here and 

below we suppress band indices) 

p q |k,a)«e J / k k+qdk '^( k ')| k + q, a ). (21) 

In other words, for small q, p q implements parallel transport described by the Berry connection A. a (k). Either 
from this observation or via a gradient expansion, we may show that at long wavelengths, the commutator 
of projected density operators at different wavevectors is 



k L 



s a (k)5>r (k+x (k_)x 7 t +c 



7k_ 



(22) 



where we define k± = k± qi + q2 . Finally, let us assume that the local Berry curvature B a (k) can be replaced 
by its average 

over the BZ; here Abz = c^/a 2 is the area of the BZ, with a the lattice spacing and Co a numerical constant 
depending on the unit cell symmetry. This yields 

[Pqi > ftfc] ~ 1 1* A P qi+q2 ■ (24) 

which is identical to the long-wavelength limit of the density algebra (20) for the LLL, with Ba^ 2 = %/2 ( ^ c '° a 
playing the role of the magnetic length £b- We may draw several inferences from Eqn. (24), as detailed in 
previous work. We summarize them briefly here. First, we observe that a coarse grained, projected, position 
operator may be defined via r cg = x cg i + y cg y = lira ^p q . It follows from Eqn. (24) that 

[X cg , Vcg] = -iB a (25) 

which identifies the r cg with the guiding center position operator in the LLL. 

11. We remind the reader that these arc simply the projected densities stripped of a Gaussian factor e~ q *b/ 4 . 
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Note that without too much work we can construct operators which obey the true GMP algebra while 
relaxing the long-wavelength approximation. For a system of N unit cells there are N points in the BZ. 
If B a (k) is truly constant we can define a set of N parallel translation operators T q for which (21) holds 
exactly: 

T q |k,a) - e^^^lk + q,^ . (26) 

The algebra of the T q is thus exactly of the form (20) without a long-wavelength restriction. We note 
that the T q are trivially isomorphic to magnetic translation operators for a system with N sites and flux 1/N 
per unit cell and it is straightforward to check that the states in the band form an N dimensional irreducible 
representation of their algebra [51]. From this perspective the idealization of a constant curvature Chern 
band hosts a algebra whose long-wavelength generators coincide with the physical density operators. 

Finally, note that both the "Hofstadter" and "Haldanc" problems give rise to (24), which unifies earlier 
lattice FQHE studies [24,25] with the ones considered here. In this sense, those earlier works can be seen 
as the first demonstration of lattice FQH phases, albeit those in which the Berry curvature is extrinsic and 
due to a magnetic field rather than an intrinsic property of the lattice model. This last observation can be 
used to get nearly constant curvature bands by approaching the Landau level limit on the lattice, i.e. by 
picking flux 1/m per plaquette and working the lowest subband at large m. While it is impossible to find a 
constant-curvature Chern band in models of Chern insulators with Af = 2 bands, it is possible to construct 
models with Af > 2 which host a Chern band with nearly constant Berry curvature [52,53]. 

Readers familiar with the LLL problem will note that there the algebra in a system with iV$ states is 
generated by N$ density operators while in a Chern band for an N x N lattice system there are N 4 states 
there are only N 2 densities (or T q if one wishes to work with a closed algebra). This distinction arises as 
the LLL is formally defined on a continuous space but is without fundamental dynamical significance as the 
relevant momenta, q£s < 1 are 0(7V$) in the LLL as well. For example, it was shown in [54] that keeping 
only this set of momenta keeps the entire physics of the quantum Hall localization transition in the LLL. 
However this counting discrepancy does have the consequence that the algebra of the densities themselves 
must close in the LLL at all q which is not the case in the Chern band. We will encounter this issue of 
operator counting again later in this review, when we discuss the Shankar-Murthy 'Hamiltonian' approach 
to the fractional Chern insulator problem. 



4.3. Geometry of the Chern Band 



The algebraic approach to the Chern band discussed can be extended further to include effects of the band 
geometry in addition to the topology [21]. This allows further constraints to be placed on models in order 
that they stabilize fractionalized phases. While important to a complete understanding of the structure of 
Chern bands, this section uses somewhat more technical machinery than the remainder, and can be skipped 
on a first reading. In the preceding section, we showed that to order q 2 , and in the uniform Berry curvature 
approximation the algebra of projected densities obeys an analog of the GMP algebra of the LLL. A natural 
extension of this result is motivated by examining higher order terms in q in the expansion of [p qi , p q2 ], and 
demanding that they in turn vanish. At order g 3 , after some algebra we can demonstrate [21] that this occurs 
if and only if a certain gauge-invariant quantity, known as the Fubini-Study (FS) metric tensor, g°j(k) is a 
constant across the Brillouin zone. The FS metric is a rank two symmetric tensor, <?g(k), which for band a 
has components [55,56,57,58] 



ft 



dh 



dkj 



dkj dkj 



E 

b 



/ f)ii a * f)n a 
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V a 6 dh 



du° 



6 dki 



(27) 



where the u"(k) are the single-particle eigenstates defined previously. It is worth pausing to expand a 
little on the meaning of this rather formidable expression. Recall that the Berry curvature of a band is 
a natural geometric quantity that emerges from considering the evolution of the eigenstates of the single- 
particle Hamiltonian h a b(k) as the parameter k evolves (adiabatically) across the BZ. Just as we defined the 
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projector onto band a we can define the orthogonal projector Q a (k) = 1 — J2p^ a 1^' &) s - > &\ ■ ^ is convenient 
to introduce the (complex) tensor 

ft£(k) = [d ki (k,a\] Q Q (k) [0 fc ,|k,a>] . (28) 

By rewriting the eigenstates in terms of the u"(k) it is readily verified that the Berry curvature B a (k) = 
— ie^Hfj (k) = — 21m [7^" y (k)] . 12 A little more work shows that the FS metric is simply the corresponding 
real part: <7g(k) = Re [7£?-(k)] . 13 Note that the gauge-invariance of both the metric and the curvature are 
manifest in the above form. 

The importance of the FS metric is that it introduces a natural notion of 'distance' in the single-particle 
Hilbert space that can induced from the k-space evolution of the single-particle eigenstates. 14 Furthermore, 
since the FS metric and the Berry curvature are the real and imaginary parts of the same complex tensor 
7Z a (k), it is also possible to show that at any k in the BZ, the trace of the FS metric is bounded from below 
by the magnitude of the Berry curvature at k, 

tr(.g«(k)) > |B a (k)|. (29) 

This geometrical constraint applies to any insulator. In the remainder of this section, we will use this in 
conjunction with the existence of a nonzero Chern number, to sharpen our understanding of the special 
properties of the restricted Hilbert space of a Chern band. 

We were led to consider the metric tensor by studying higher order terms in the expansion of the algebra 
of projected density, and we identified its uniformity as a criterion for this to match the algebra of 
the LLL to order q 3 . Thus, we have found an additional criterion for identifying "good" band structures 
from the point of view of hosting interacting topological phases. A natural question to ask is whether one 
obtains an infinite set of such constraints, order-by-order in q; it would be somewhat deflating if this were 
the case. Remarkably, we will demonstrate that when the band structure satisfies one additional constraint, 
the Chern band projected densities satisfy the algebra of projected LLL densities at all orders in q. 

To this end, it is useful to consider the transformation properties of TZfj under unimodular transformations 
of the coordinates, i.e. those of the form x a — > x' a = A a bXb with det A = 1. The FS metric transforms in the 
usual way, (<7"J'(k) = A ac A(,d<7" d (k), while the Berry curvature is invariant. Moving to 'primed' coordinates 
in which g Q (k) is diagonal, 15 using the invariance of the determinant under unitary transformations and 
applying the inequality (29) we obtain 

2v/dct 3 «(k) = 2 A /dct(< 7 °(k))' = tr($(k))' > |[S Q (k)]'| = \B a (k)\ (30) 
Since this is true across the BZ, we can integrate this to find 

/ d 2 k dot g a (k)>] [ d 2 k\B a (k)\ 2 >]A BZ B 2 a = ^ (31) 

JBZ 4 JBZ 4 A BZ 

where we use the fact that the square-averaged Berry curvature is bounded from below by the square of the 
Chern number per unit area of the BZ. Eq. (31) is the key relation between the Chern band geometry to 
its nontrivial topology: namely, the integral of the determinant of the Fubini-Study metric is bounded from 
below by a number proportional to the square of the Chern number of the band. 

We now consider the case when (i) the inequality (31) is saturated and (ii) the FS metric is uniform in 
the BZ. It is easily seen that these conditions amount to requiring, in addition to the uniformity of the 
Berry curvature obtained previously, the single additional constraint that det(g"(k)) = ||,S Q (k)| 2 at all k 



12. Note that the curvature also has indices, but in d = 2 there is only one nonzero component so we can suppress them. 

13. Observe that owing to the equivalence between the position operator Xi and the k-spacc derivative —id/dki it is possible 
to rewrite = (k, a\xiQ(k)xj |k, a); in this form we see that the non-vanishing of the Berry curvature roughly corresponds 
to the non-commutativity of the x- and y-componcnts of the position operator. 

14. This is intimately connected to the geometry of Kahler manifolds as well as to the geometric theory of the insulating 
state, but unfortunately we cannot discuss these beautiful subjects here [59]. 

15. For d = 2 in this basis simple algebra shows that the square of the trace is four times the determinant. 



16 



in the BZ. In this limit and in the 'primed' coordinate system used previously, in which the FS metric is 
diagonal, we have [<7^(k)]' = ^f-Sij. Writing the density operators in the rotated coordinate system, it is 
straightforward to demonstrate that they satisfy a generalized, metric-dependent version of the Woo algebra: 

[fta,*,] = 2*sin (SI^El) e^^p^ = 2*sin e*^' ^ (32) 

There are three important aspects of this result worth reflecting on. First, note that Berry curvature and 
the FS metric both appear in the form of the Woo algebra obtained in (32), which thus algebra also applies 
to bands which have a higher Chern number and therefore differ fundamentally from Landau levels with 
C = 1. Second, observe that the conditions under which we get a closed algebra of the projected density 
operators can be stated purely in terms of the FS metric. Third, we observe that for the C = 1 case, 
(32) looks remarkably similar to the algebra of projected densities that obtains in a Landau level, with the 
combination e^ l9lm ^ m playing a role analogous to the LL-dependent form factor that relates projected 
densities to the magnetic translation operators (which are identical in every Landau level) for the continuum 
FQHE. A heuristic interpretation of the additional constraints is therefore that they quantify not simply how 
close Chern band projected densities are to realizing the continuum Landau level algebra, but in addition 
how close they are to realizing the lowest Landau level. 

For a system where the ideal conditions under which this algebra is obtained do not hold, the degree of 
deviation from these conditions provides an additional criterion to quantify how favorable a Chern band 
is for hosting FQH-like physics. Conversely, if one finds fractional topological phases in systems where the 
deviations from these conditions is considerable, one could argue that the physics of those systems is new 
and distinct from the conventional fractional quantum Hall effect. The effects of disorder also enter the 
Hamiltonian through terms that involve the projected density operator. This suggests that the effects of 
disorder in the Chern band are likely to be the same as in the LLL when the conditions stated above for 
the FS metric are satisfied. Extensions of the algebraic approach to higher dimensions and to time-reversal 
symmetric Z 2 topological insulators have also been discussed in the literature [60,61,62]. 

4.4. Operator Counting, Smoothed Densities, and the Shankar-Murthy Approach 

We now return to the question of operator counting that we discussed earlier. As motivation, we first 
observe that in cases where the deviation of the Berry curvature from its average value is bounded, |B a (k) — 
B<x\ < \Ba \ ~ e ; we mav define a 'smoothed' density operator which may be regarded as the projection of an 
operator p s (r) local in position space; for ga< 1 this gives a modified form of (21) 

^ M = ^kT e4/kk+qdk ''" 4a(k ' )|k+q ' a) - (33) 

At long wavelengths, the algebra of smoothed densities closes, and in this limit (24) is an exact equality 
when p q is replaced by p*. This provides a clue that there are observables in the Chern band that indeed 
satisfy a well-behaved CMP algebra, related to the projected density operators in some nontrivial fashion. 

Building on the different operator counting in the LLL versus the Chern band, Shankar and Murthy 
were able to take a slightly different perspective on the operator algebra approach to FCIs. Their approach 
(described in Ref. [23,22]) consists of three separate steps: 

(i) They construct a set of one-body operators PGMp(q + G) which act only within a Chern band, and 
satisfy the full GMP algebra (20). Here, q is a crystal momentum (restricted to lie in the first BZ) while 
G is a reciprocal lattice vector. For an TV x N lattice, the number of independent crystal momenta is N 2 , 
and Shankar and Murthy demonstrated that the 'umklapp' operators PGMp(q + G) are independent 
only for N 2 choices of reciprocal lattice vector G. Thus, together there are N 4 independent GMP 
operators PGMp(q + G) which thus form a complete basis in terms of which one may express any of 
the N 4 one-body operators in a Chern band. 

(ii) They then demonstrate that the projected density operators in the Chern band may be rewritten in 
terms of the PGMp(q + G): 
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Figure 4. Mappings within and between the LLL and the Chern band. Path 1: Within the LLL, there exist special 
'model Hamiltonians' designed to a render a certain FQH state an exact ground state; a wealth of numerical evidence suggests 
that in many cases these are adiabatically connected to the ground states of more realistic Hamiltonians, for instance those 
obtained by projecting a density-density interaction to the LLL. Path 2: Qi and Shankar-Murthy both construct mappings 
from the LLL to the Chern band that allow them to construct model Hamiltonians in the Chern band, with the same topological 
order in the LLL. Path 3: The present authors instead considered realistic density-density interactions projected to the Chern 
band, and gave conditions for this to be adiabatically connected to the FQH ground state of a similarly realistic Hamiltonian in 
the LLL . Path 4: Recent numerical studies suggest that the model Hamiltonians are adiabatically connected to these more 
realistic projected density-type interactions. 



where the expansion is over reciprocal lattice vectors and the coefficients are determined from data 
on the lattice realization of the original Chern band. Thus the projected densities are a specific linear 
combination of operators satisfying the GMP algebra - which explains why the /? q themselves only 
satisfy (20) in an idealized limit. They also noted a distinction between the projected density and the 
exponential of the guiding center density in a Chern band (in the LLL these are the same, up to a 
normalization depending only on q). Recall that an analogous situation arose in our discussion of the 
Chern band geometry in the preceding section. 

(iii) The preceding observations allow them to express a Hamiltonian projected to the Chern band in terms 
of operators that obey the GMP algebra. Note that since one may naturally formulate the full LLL 
problem in terms of the GMP operators, this suggests another route to constructing model Hamilto- 
nians for a given FQH phase. Furthermore, Shankar and Murthy demonstrate that the 'Hamiltonian 
theory' that they pioneered for the FQH - in which Chern-Simons flux attachment is placed on an 
operator footing - can be ported to the Chern band by working with the pgmp (l) • This in turn permits 
the computation of gaps and response functions within a Hartree-Fock mean field theory. 

For reasons of brevity we cannot give a detailed elucidation of these points, and encourage the reader to 
study the original articles [23,22] for a pedagogical discussion of the Hamiltonian theory of fractionally filled 
Chern bands. 

4.5. Adiabatic Continuity and Establishing Topological Order 

We have provided three different analytical perspectives on the correspondence between Chern bands and 
the (lowest) Landau level. We now place these in context by describing how they are related to the problem 
of establishing the topological order of the ground state of an interacting, fractionally filled Chern band. For 
this, it is useful to recall that at a filling for which a quantum Hall state is known to exist, there are typically 
two points in the space of all LLL-projected Hamiltonians that are of especial interest. The first is the model 
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interaction, which we denote -ff^odoi' that ren ders a certain FQH wavefunction an exact gapped ground 
state 16 . The second is the one that is reached by projecting a realistic electronic interaction Hamiltonian 
into the LLL, which we shall call -f^^aiistic- 17 There is much numerical evidence establishing that in many 
cases the ground states of H^^ istic are connected adiabatically to those of i?modeii tnus establishing that for 
experimentally relevant and reasonable interactions the ground state has the same topological order as that 
encoded by the model system. 

What about the Chern band? Here, we must first produce model Hamiltonians H^ dcl on the lattice, since 
those constructed in the continuum rely intimately on the analytic structure of the LLL wave functions and 
cannot be ported directly to a Chern band. Here, Qi's prescription for mapping Landau gauge eigenstates 
from the LLL to the Chern band provides one route. Another is that of Shankar and Murthy; by suitably 
picking a pseudopotential interaction in the LLL and replacing the corresponding projected density operators 
by the Pgmps, they can construct model Hamiltonians with similar topological properties in the Chern 
band. However, both these approaches can lead to Hamiltonians which may look somewhat unnatural on 
the lattice - Qi's construction because the Wannier states may not be related by a simple transformation to 
the densities, 18 they are not and Shankar and Murthy's approach because of the nontrivial relation between 
the pgmpS and the projected density. 

A different approach, suggested by the present authors, attempts to simply study the properties of 
-^realistic' defined by projecting a density-density interaction of some fixed range to the Chern band. Here, 
the mismatch between the CB and the LLL is encoded in the structure of the algebra of projected den- 
sities: as we have emphasized, they obey the GMP algebra only in certain idealized limits. Nevertheless, 
this suggests that when conditions are such that this idealization is not too far-fetched, one can argue that 
the ground states of ff^aiistic anc ^ Realistic wm be- in the same phase, and draw on the continuity between 
-^realistic an d ^modoi m the continuum Landau problem to establish the topological order of the Chern band 
ground state for realistic interactions. Alternatively, one may view the idealization as placing constraints on 
various parameters, such as Berry curvature fluctuations, for there to exist an adiabatic path in the Chern 
band connecting -Holistic to H^ dcl . We note that numerical studies have demonstrated that when these 
constraints are violated, the stability of the Chern band ground state is adversely affected; also, the adiabatic 
continuity between Qi's model Hamiltonians, and more realistic interactions, has recently been verified [48]. 

5. Parton Constructions and Phase Transitions 

Thus far, we have focused on matching the physics of a Chern band with that of the lowest Landau 
level, and in establishing the existence of lattice analogs of familiar FQH states. Our survey of the existing 
analytical and numerical evidence should at this point have convinced the reader that at least some of the 
phases seen in the lowest Landau level survive to the lattice limit. In this section and the next, we change 
perspective slightly and ask whether the presence of the lattice can lead to additional physics beyond that 
familiar from the Landau level. A natural method to analyze such question is to use the 'parton' construction 
or projective construction, a mean field theory that has had considerable success in analyzing both continuum 
FQH phases as well as lattice spin liquid models. As has been our theme overall, we will attempt only to 
give a flavor of this vast area, and identify a few interesting results that have emerged from such studies. 

5.1. Partons and Projective Symmetry Groups 

We will give a very simplified account of a parton construction [63] of the v = 1/m Laughlin state, first in 
the continuum and then its lattice analog. In the continuum, the basic observation is that as a consequence 
of the dependence of the quantum of flux on the electric charge, a Landau level that is at 1/m filling for 



16. Note that not all FQH states have parent Hamiltonians of this form; here we shall restrict ourselves to those that do. 

17. Of course, there may be other details, such as the range of the interaction, etc., but we shall ignore this subtlety and 
discuss all 'realistic' Hamiltonians in a unified fashion. 

18. Recall that in a Chern band, Wannier orbitals can be chosen to be eigenstates of only one position coordinate. 



19 



particles of charge e, is fully filled for particles of charge e/m. In other words, an identical number of e/m- 
charged particles would have a unique, Slater determinant ground state. In order to put this observation to 
use, we divide each electron into m fcrmionic 'partons', each of charge e/m, distinguished by a flavor index, 
i.e. we write the electron operator as 

m 

c(r) = J] fa(r) (35) 

Q = l 

We can now write mean-field states for each parton flavor independently. However, in splitting apart the 
electron and treating the resulting partons independently, we artificially enlarged the Hilbert space. To 
amend this, the final step of the construction is to project the mean-field wavefunction back to the physical 
subspace, by enforcing the requirement that any valid many body wavefunction for the physical electrons 
have m partons at the same position coordinate. In other words, we write for the electron wavefunction 

N m 

*y ro ({rj)=(oinn^wi MF )- ( 3g ) 

i— 1 a— 1 

Implementing this procedure for the LLL in the disc geometry, and using the fact that the charge e/3 partons 
fill Landau levels as discussed above, we find that 

*V"* ({r,}) = [^=i on ( { rj)] m = H( Zl - zjre-Xt (37) 

where z = x + iy, £b = \Jhje\B\ is the electron magnetic length, and we used the fact that a filled LLL 
of charge e/3 partons has the Slater determinant wavefunction parton ({ r i}) = Ili<j( z i — z j) e ~ ^ B > 

where is — \/ml B is the magnetic length of the partons. It is useful also to consider a field-theoretic 
formulation of the parton construction. The LLL parton decomposition (35) of the electron has a SU(m) 
gauge symmetry: rotations in the internal space that mixes the partons yield the same electron operator, and 
such rotations can be performed independently for each r. Essentially, within the field theory the fluctuating 
gauge field emerges as a way of enforcing the local constraint that glues together partons to form physical 
electrons. The emergent gauge symmetry is reflected by the symmetries of the mean-field Hamiltonian that 
describes the parton, which is simply three copies of the kinetic energy of a charge e/m particle in an external 
field; since the partons are gapped (fill a Landau level), they may be integrated out, yielding an SU{m)\ 
Chern-Simons theory, that yields the correct ground-state degeneracy on the torus, and other topological 
properties. Note that the assumption of full SU{m) symmetry is not strictly necessary; it is possible also 
to consider terms in the mean-field Hamiltonian that mix partons and break the symmetry down to the 
discrete Z m subgroup, but in the restricted Hilbert space of the LLL, both these symmetries give rise to the 
same electronic wavefunction after projection. 

How does the parton approach fare on the lattice? Our discussion will follow that of Lu and Ran [63] , who 
first applied the parton construction to Chern insulators, but once again we only give the broad outlines. 
We begin by performing a similar decomposition (35) of the electron operator on each lattice site into three 
partons, assuming for simplicity that each site has only one orbital (but of course, multiple sites within each 
unit cell.) A generic mean-field Hamiltonians for the partons takes the form 

m 

h mf = ^2Y1 fl(r)M a0 (r,r')Mr) (38) 

r.r' a. [3—1 

where M(r, r') is an m x m matrix. Under a local SU{m) gauge rotation that transforms partons via / a ( r ) — > 
Gap( r )fp( r ) M transforms as M(r,r') — > G(r)M(r, r')G(r'), with multiplication in the parton indices 
implied in the latter expression. Once again, the final electronic wavefunction is obtained by diagonalizing 
the quadratic parton Hamiltonian (38) and projecting to the physical degrees of freedom by imposing the 
constraint that m partons must always appear togthcr, once again using (36). We turn now to the delicate 
question of how to choose H MF to describe a fractional Chern insulator. 
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As motivation, it is useful to note that for a v — 1/m Chern band, there are is on average 1 electron per 
m unit cells; thus, a mean-field electronic wavefunction will generically lead to a metallic state. Ideally we 
would like each of the to partons to fill a band with Chern number 1; this would replicate our construction 
in the continuum. However, in the latter case, there was a simple expedient, charge fractionalization, which 
allowed us to change the effective filling via a redefinition of the quantum of flux for the partons. Here, 
the presence of the lattice makes this more subtle: essentially, the increase in the quantum of flux for the 
partons implies that, although the electron band structure respects all the lattice symmetries, the partons 
have a fraction of a flux in each plaquette of the original lattice. How are we to construct a symmetric 
wavefunction? 

The key lies in the notion of the projective symmetry group familiar from slave-particle theories of spin 
liquids [64]. In essence, the physical degrees of freedom are the electrons, so that a broken symmetry for 
the partons that is invisible to observables written in terms of electron operators is irrelevant. For instance, 
two mean-field solutions that differ by gauge transformations are equivalent as they yield identical electronic 
wave functions after projection; similarly, mean-field states may break lattice symmetries, and yet give rise to 
symmetric electronic wave functions. In mathematical terms, the various mean-field ansatze form projective 
representations of the symmetry group, and can be classified using well-developed technology. 

Following this general line of reasoning, the path to constructing symmetric wave functions at v — 1/m is 
as follows: we obtain a mean-field parton band structure by enlarging the unit cell on the original lattice, for 
instance by sticking 2ir/m flux through each original unit cell. Since there was originally one electron per to 
unit cells of the original lattice, but in the mean-field description each electron has been decomposed into 
to partons, and the new unit cell consists of to of the original cells, it follows that the mean-field state we 
construct in this fashion has to partons filling the m lowest bands, and the resulting parton state is gapped. 
Note that counting in terms of the electrons, this corresponds to one electron per (enlarged) mean-field 
unit cell so that we expect the gap to be robust under the projection back to electronic degrees of freedom. 
Since each parton has charge e* = e/m, and there are to filled C = 1 bands, we obtain the desired Hall 
conductance 

(e*) 2 1 e 2 

o- xy = m— — = — — - (39) 
h m h 

Although we have assumed each of the to parton bands has C = 1, note that different choices can give states 
where the filling and the Hall conductance do not bear a simple relation to each other [39,23,22]. It is also 
instructive to view the construction in terms of an explicit wavefunction. The electronic wavefunction is the 
projection of the Slater determinant for to filled parton bands: 

*e£l ({*» = ^det W; W {a - 1)N+i , {n - 1)N+j = v%\rf) (40) 

where V implements the constraint, and W is an mN x mN matrix, with i,j — 1,2, ... , N and a, n = 
1,2,..., to. Here, u^irf), is a single-particle momentum-kj state in which parton a is in the n th band 
(counting from below) obtained by diagonalizing the mean-field parton band structure. When the mean- 
field band structure has full SU{m) symmetry, M a p(r,r') — M(r, r')<5 Q ^, the lowest to bands of the parton 

band structure are identical and we have u^(rf) — (r^)5 n!(x . The general form (40) then factors into 
the rn-th power of a Slater determinant: 

K^u {m) = [det[u^i)]] m - (41) 

The fully SU(m) symmetric case is thus analogous to the case of the LLL. However, in general, tunneling 
terms between the partons are present and rule out such a factoring, restricting us to Z m parton wave 
functions of the form (40). This distinction between the Z m and SU{m) wave functions is one way in which 
the parton construction on the lattice differs from that in the continuum LLL. 

From the preceding discussion it should be clear that there is quite a bit of freedom in implementing this 
procedure; ideally we would like to be able to exhaustively list all such (mean-field) wave functions. Here, 
the full technology of projective symmetry group (PSG) classification is required, and we refer the reader 
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to [63], as well as related work [65,63,27] for details. We content ourselves with noting that the presence of 
the lattice allows the possibility of 'symmetry-enriched' topological phases, whose low-energy descriptions 
have similar topological structure but realize lattice symmetry in distinct ways. 

5.2. Bandwidth- Tuned Phase Transitions 

A fractionally filled Landau level in the continuum limit has no unique ground state in the absence of 
interactions. In contrast, a Chern band (alternatively, a Landau level with a strong periodic potential) 
typically has a nonzero dispersion, characterized by its bandwidth t (while perfect flat bands have t = 0, 
these are fine-tuned.) This weak dispersion can select a ground state in the U/t — > limit, where U is as 
usual the interaction scale; at partial filling, this ground state is expected to be compressible. 19 As U/t is 
increased, eventually we enter the flat band limit where the kinetic energy can be ignored; we assume the 
system forms an incompressible FQH ground state in this limit. A critical point separating the compressible 
and incompressible phases is a possibility unrealized in the ideal limit of the FQHE. Although such a Mott 
transition in a Chern band could always be first-order, there is the additional intriguing possibility of a 
continuous transition into FQH states. 

Several examples of such transitions at half-filling were studied by Barkeshli and McGreevy [27,28] and 
can be divided into three sets of phases which can have continuous transitions between each other: 

(i) between superfluid, Mott insulating and v = 1/2 Laughlin states of bosons; 

(ii) between an ordinary Fermi liquid, a gapless Mott insulator, and a composite fermion Fermi liquid (i.e., 
the Chern band analog of the Halperin- Lee- Read state in the half-filled Landau level); and 

(iii) between a p + ip paired fermion superfluid, the B-phase of Kitaev's honeycomb model, and a lattice 
analog of the Moore Read Pfaffian phase. 

All these were studied within a parton mean-field theory which is expected to capture the gross features of 
the phases as well as the critical behavior between them. Their results may be summarized as follows. In 
essence, through the parton decomposition, transitions between the different phases can be rewritten in terms 
of Chern number- changing transitions between insulators: | AC| = 1 between successive members of each of 
the three sets above. A general critical theory for Chcrn-number changing transitions in two dimensions can 
be formulated in terms of massless Dirac fermions; a change in the sign of a Dirac mass changes the Chern 
number by 1. In the cases studied by Barkeshli and McGreevy, the massless Dirac fermions are also coupled 
to a Chern Simons gauge field, which emerges from the parton construction. For all three sets of phases, 
a direct transition between the first and the third phase requires |AC| = 2 and is generically multicritical, 
unless an additional symmetry, such as inversion, is present to force two Dirac fermions to undergo a 
simultaneous change in the sign of their masses. 20 Furthermore, the parton construction describing (i) is 
the parent example; the remaining two follow from a further parton decomposition in which an additional 
neutral fermionic degree of freedom is attached to each boson in (i) and allowed to form either a Fermi 
surface of gapless excitations (ii) or a paired state (iii). We direct the interested reader to Refs. [28,27] for a 
detailed discussion of these transitions, along with a list of physical manifestations in crossover physics and 
scaling behavior. 

6. Higher Chern Numbers and Nonabeliana 

Another possibility that can be realized in lattice models is a band with higher Chern number; the interest 
in these examples is primarily concerned with the possibility of realizing non-Abelian statistics in a manner 
intrinsically tied to the lattice limit of the FQHE in a Chern band. In this context we note two intriguing 
directions we believe are particularly worthy of further attention. The first is a numerical study by Zhang 



19. Note that this is the expectation whether the particles are bosons (a superfluid state) or fermions (a filled Fermi sea). 

20. An alternative possibility that gives AC = 2, a quadratic band-touching, is perturbatively unstable in d = 2 and thus 
ruled out. 
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and Vishwanath [29] , who analyzed wave functions obtained by taking N copies of a Slater determinant for 
a fully filled Chern number C band and Gutzwiller projecting to a Hilbcrt space that obeys a single-site- 
occupancy constraint. A careful variational Monte Carlo study of the entanglement entropy of degenerate 
ground states on a torus for the case N = C = 2 provides strong evidence that the resulting wavefunction 
has the same topological order as the Moore-Read Pfaffian state, described by an SU(2) 2 Chern-Simons 
theory. Combined with a previous analysis of lattice analogs of bosonic Laughlin states along similar lines 
[66], this suggests that, in general, projecting the N th power of a Slater determinant of a filled Chern number 
C = k band provides a lattice version of a quantum Hall state whose topological order is identical to that of 
the SU(N)k Chern-Simons theory coupled to fermions (which also describes a system of k 'layers' at total 
filling v = k/N.) 

The interplay between the lattice-scale physics and topological order in a band with C > 1 was studied 
in detail by Barkeshli and Qi [67]. Their result can be readily understood within the Wannier description; 
recall that the center-of-mass position of a Wannier function shifts from x — > x + C when k y evolves from 
to 2-7T. If C > 1, it is therefore possible to define C distinct families of Wannier functions that evolve 
independently under k y — > k y + 2ir. For concreteness, consider the C = 2 case. Here, the Wannier functions 
on even and odd sites form distinct families: in the notation introduced previously, we define 

|W^ s=fcB+2wn ) = |W(fc W) 2n-l)), \W 2 Ky=ky+2nn ) = \W(k y ,2n)) (42) 

On its own, the observation of two families of cigenfunctions is not particularly remarkable; a similar structure 
obtains, for instance, in multicomponent quantum Hall systems such as quantum Hall ferromagnets and 
bilayer systems. The key difference is the action of the lattice symmetry operations on the 'internal' index. 
For instance, the translations along x interchange the internal index, which are invariant under translations 
along y: 21 

f x \W 1 Ky ) = \W 2 Ky ), f x \W 2 Ky ) = \W 1 Ky ), f y \W a Ky )=e iK v\W a Ky ),a=l,2 (43) 

As a consequence, lattice defects - such as dislocations - have a nontrivial impact on the underlying topolog- 
ical state. In a concrete example, Barkeshli and Qi [67] considered Abelian 'bilayer' FQHE states realized in 
a C = 2 band, and showed that such defects that act nontrivially on the this internal 'Z2' degree of freedom 
have a nontrivial topological degeneracy - analogous to quasiparticle states in a non-Abelian state. In effect, 
by acting nontrivially on the layer index, these 'twist' defects simulate the action of a change of genus of 
the underlying space - and are hence dubbed 'genons'. A subsequent detailed analysis [69] of the projective 
non-trivial braiding statistics obtained in this fashion, while beyond the scope of the present article, suggests 
an interesting new line of enquiry with possible implications for topological quantum computation. 

7. Proposed Realizations 

Thus far our discussion has focused on lattice models and numerical results. As we emphasized earlier, 
an appealing feature of the lattice FQHE phases discussed above is that the characteristic energy scale 
at which new physics sets in is tied to the microscopic (lattice) scale, which is quite large; if they can 
be experimentally realized, it is possible that Chern bands will exhibit fractionalization and associated 
phenomena at significantly higher temperature scales and be considerably more robust than in 2DEGs. We 
flag a few especially promising examples below. 

In a solid-state context, Xiao and collaborators suggested that topological bands may be engineered 
in oxide interaces [30] . From a combination of first- principles band structure calculations and tight-binding 
modeling, they predict that bilayer [111] heterostructures of transition-metal perovskite oxide host relatively 
flat bands with a nontrivial Z 2 topological index. Apart from providing a route to room-temperature quantum 



21. Recent work [68] has shown that a slightly more involved definition of the Wannier orbitals and the associated action 
of the translation operator is required in order to impose boundary conditions that fully respect various lattice symmetries. 
Nevertheless, the interplay of lattice defects with topology and the field theory describing this structure remain, so we will 
ignore these subtleties. 
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spin Hall insulators, when the bands are split by a Zeeman term (this can be extrinsic, or in an interacting 
model emerge due to ferromagnetic ordering) they emerge with a non-zero Chern number. The flatness, 
which is due to the weak dispersion of the crystal-field split rf-orbitals, suggests that fractional quantum 
Hall physics is a plausible scenario in these systems; see also [36]. 

Another direction is to directly engineer spin-orbit coupled models in optical lattices of ultracold atoms. 
Cooper [31] has shown that the strong-lattice limit can be achieved in two-species ultracold atom systems 
(where for instance the two species are the ground and long-lived excited states of an alkaline earth atom) 
effectively described by the single-particle Hamiltonian 

l+V s (r)-<r (44) 

where the 'spin' er refers to the level index. Here, Vo(r) is a conventional scalar potential term, for instance 
from a conventional optical lattice, while V s (r) reflects a laser pattern that generates an 'optical flux 
lattice', where the diagonal components reflect a species-dependent scalar potential, and the off-diagonal 
terms represent a Raman coupling that tunes the rate of interspecies conversion as a function of position. 
A judicious choice of laser parameters (Vb, V s ) leads to a lowest band with a non-zero Chern number. It 
is also possible to adjust parameters such that the band is relatively flat [70]. This and similar proposals 
suggest that the high range of tunability afforded by optical lattice models will allow us to study a variety 
of features of Chern bands. 

A different bosonic cold atom implementation is to consider the rotational degrees of freedom of ultracold 
molecules, whose actual spatial motion is pinned by a deep optical lattice. Specifically, Yao et. al. [32,33] 
consider a system of such pinned, three-level dipolcs in a two-dimensional lattices. A fractional Chern 
insulator Hamiltonian can be shown to describe the spin flips of this problem which are effective hard 
core bosons. The time reversal breaking necessary for a nonzero Chern number is provided by elliptically 
polarized, spatially modulated electromagnetic radiation; the nontrivial spatial dependence of the hopping 
relics on the the long-range dipole-dipole interactions characteristic to polar molecules, which naturally 
introduces orientation-dependent hopping amplitudes. Appropriately tuning these also allows control of the 
hopping interference and permits bands with / > 10 [33]. Numerical studies of these models show a variety 
of phases at v = 1/2, including fractional Chern insulators for a range of parameters that arc not too far 
from current experimental capabilities [32]. 

8. Concluding Remarks 

We hope that the reader who has persevered to this point has acquired an impressionistic perspective 
of several quite different aspects of this rather active area of research. We have taken a more or less his- 
torical approach, first explaining how to design topological flat band models and providing an overview 
of the numerical results, before drawing analytical parallels between fractional Chern insulators and FQH 
physics in the lowest Landau level. We then briefly touched upon a host of interesting topics: construct- 
ing and classifying fractional Chern insulator wave functions using parton mean-field theories, the theory 
of bandwidth-tuned transitions into FQH phases, and various aspects of lattice physics in bands of higher 
Chern number. Each of these has spawned its own numerical and analytical offshoots, and we have attempted 
to provide pointers to the relevant literature wherever possible. We also gave an overview of recent progress 
in developing experimentally realistic scenarios where these correlated phases emerge. 

In closing this review, we take a somewhat longer view, and outline several directions which seem (to us, 
at least) worthy of further study, although we hope that the preceding discussion has already prompted our 
readers to form similar questions independent of these. Once again, we sound a note of caution: given the 
rapid progress in the field, it is quite likely that this list of future directions might quite rapidly become 
dated as new results emerge on various fronts. 

Disorder. — The first of these is the important and (at this point) relatively little-studied problem of 
disorder in the Chern band. Here it is worth noting that in addition to the square root of the mean Berry 
curvature (which as we have seen plays a role analogous to the magnetic length in the LLL), we have 



H = 



2m 
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a second (inverse) length scale: the characteristic momentum-space scale k a for variations of the Berry 
curvature. Naively we might expect that the role of disorder depends on how its strength - as quantified by 
the mean free path, ^mf ~~ compares to one or both of these scales. For instance, in a problem with sufficient 
disorder, most notably the problem of localization within the CB, we expect that impurity scattering will 
effectively lead to an averaging of the curvature over the band and give results similar to those in the LLL. 
This is consistent with what is known about the (integer) QH transition in Chern insulators [71] although a 
direct test in the projected CB would be desirable. A simple estimate for the requisite strength of disorder 
can be made by comparing Imf to fc" 1 ; for disorder sufficiently strong that Imf ~ k~ 1 , the random potential 
will scatter between points in the BZ that are k a apart, and thereby average their Berry curvature. More 
ambitiously, it would be quite instructive to construct a global phase diagram for fractional Chern insulators, 
akin to that proposed by Kivelson, Lee and Zhang [72] for the FQHE in the LLL. Here specifically it is worth 
recalling that the full global phase diagram for the continuum Landau problem emerges only when mixing 
between different LLs is taken into account; for fractional Chern bands it is possible that the mixing between 
Bloch bands leads to qualitatively different physics. 

Long-Range Interactions and Competing Phases. — Turning for a moment to the continuum problem, we 
observe that there fractionalizcd phases emerge in a comparatively narrow window of parameter space: they 
compete with integer QH liquids as well as electronic crystalline phases, the latter in particular drawing 
enhanced stability from long-range Coulomb physics. We note that so far numerical studies in the Chern 
band have been largely confined to Hamiltonians with relatively short-ranged interactions. For the FQHE 
this is typically justified by invoking screening to argue that the Coulomb repulsion is operational over no 
more than a few magnetic lengths - although we note that even in the LLL there are well-known examples 
where the long-range tail of the Coulomb interactions carries a sting [73]. In a CB in contrast, the only 
length scale is the inter-atomic spacing, and even a screened interaction can be quite long-ranged on the 
lattice scale. The role of such longer-ranged interactions in stabilizing or inhibiting the formation of FQH 
phases, and in particular in supporting competing crystalline phases, may be quite relevant to experimental 
systems. We note that some numerical work has begun to explore competition between phases in the Chern 
band [74]. 

Edge States. — An integral aspect of the theory of the FQHE is the Luttinger liquid description of the 
gapless edge modes [75]. The structure of the edge of a fractional Chern insulator has received comparably 
little attention. While presumably the universal edge properties are captured by the standard chiral Luttinger 
liquid theory, it is possible that the ability to produce a lattice-scale edge and the interplay of lattice 
symmetries with the edge may provide a source of additional complexity, worth exploring. 

Experimental Detection. — As we have described already, experiments on fractional Chern insulators 
remain in their infancy, although several proposals seem quite compelling and suggest that a realization even 
with existing technology is not too far off. In this context it is important to develop techniques to characterize 
the microscopic properties of topological band structure, which as we have seen play an important role in 
stabilizing correlated phases. In light of this, a recent proposal [76] to measure the Zak phase in Bloch 
bands using a combination of Ramsey intereference and Bloch oscillations are promising, and have been 
demonstrated in experiments on ID optical lattices [77] . Equally important is the development of methods 
to identify correlated topological phases, particularly since the presence of the lattice leads to ambiguities 
in interpreting standard measures such as the Hall conductance. 

Matrix Product and Tensor Network Representations. — The study of gapped phases in one dimension 
has benefited greatly from the observation that they have ground states that can be efficiently represented 
as 'matrix product states' [78]. The power of the matrix product form is twofold: first, the density-matrix 
renormalization group (DMRG) algorithm [79,80] efficiently obtains the best matrix product representation 
of a ground state of a given local Hamiltonian; and second, the structure of the MPS allows efficient analysis 
of its behavior under symmetry operations. Similar tensor network constructions exist in higher dimensions 
[81] but a full understanding of these is still being built. The FQH phases are gapped correlated phases, 
and recently several groups have reformulated the model wavefunctions as matrix product states [82,83]. It 
would be quite instructive to develop a similar understanding for FQH phases in Chern bands, as well as to 
develop efficient algorithms to simulate them. 

Extensions. — We turn finally to the important and interesting question of what other correlated topolog- 
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ical phases might emerge in lattice models. One obvious extension that has already received quite some nu- 
merical and analytical attention is the problem of constructing interacting analogs of time-reversal invariant 
topological insulators in two and three dimensions. Unlike the FQHE, this is a problem with little or no prior 
history, so there are many open questions; some work along these directions already exists [84,85,86,87,88,89]. 
A more recent development is the classification of 'symmetry-protected topological phases' (SPTs) [90]. These 
are gapped phases which have no topological order (they lack ground-state degeneracy and fractionalizcd 
excitations) but still host gapless edge states, and cannot be adiabatically continued to a trivial gapped 
phase as long as a certain symmetry is left unbroken - for instance, topological band insulators of non- 
interacting fcrmions may be thought of as SPTs protected by time-reversal symmetry, while the Haldanc 
phase of a spin-1 chain is protected by 7r-spin rotation, time reversal, or inversion symmetry[91,92] . 22 Since 
the understanding of SPTs grew out of very similar considerations as the integer and fractional QHE, it 
is quite natural to consider Chern bands as a natural venue where they might be realized. A final note of 
complexity emerges when one considers the interplay of lattice symmetries and topological order, as in a 
few examples above; work in this direction has potential to reveal signatures of hard-to-detect topological 
phenomena in more conventional probes of the solid state. 
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